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We study pion curvature radiation by a proton, i.e. pion emission by a proton moving along 
a curved trajectory. We suggest an approximate semiclassical solution and the exact solution for 
which we assume that a proton moves in a fictitious magnetic field with the Larmor radius equal to 
the curvature radius of the real trajectory. As possible application we consider the pion radiation 
by ultrahigh energy protons moving along curved magnetic field lines. Such situation can occur in 
the magnetosphere of a young pulsar, in the magnetosphere of the accretion disk around a black 
hole, and in the vicinity of a superconducting cosmic string. The decay products of these pions, 
such as high energy photons or neutrinos, can give the observable consequences of the considered 
mechanism. 



I. INTRODUCTION 



According to the commonly used terminology, curvature radiation is electromagnetic radiation by a charged particle 
moving along a curved trajectory. Naturally, the implied mechanism can be generalized to the emission of any kind 
of particles. A simple example is pion curvature radiation by a proton. Protons can move along a curved trajectory 
in the gravitational field, e.g. in the vicinity of a black hole. Another example is the movement of a proton along a 
magnetic field line in the case of a strong field, e.g. in the magnetosphere of a pulsar or of an accretion disc. The 
transverse momentum of a proton in this case is suppressed due to synchrotron radiation. 

This formulation of the problem is semiclassical by its nature: we define the interaction between particles in a 
quantum way, but introduce at the same time the classical concept of a trajectory, which does not exist in quantum 
field theory. We suggest here some solutions to the problem for the case of 7r° emission by a proton and indicate 
possible implications. 



II. KINEMATICS AND DYNAMICAL PARAMETER OF THE PROBLEM 



To make a quantitative insight into the problem let us consider a proton moving along a curved trajectory with 
local curvature radius Rc- We introduce the coherent interaction length, knt, along the trajectory on which the pion 
is radiated. The four-momentum of the proton at the begining and the end of Imt is denoted by {Ei,pi) and (Ef,pf) 
respectively, and that of a pion by (e,fc). We describe the effect of the curved trajectory in the interaction region 
by the longitudinal and transverse (relative to pi) momenta and q± given by some external field and neglect the 
change of proton energy (\q\ ^ \p\). Then we have 

Ef ^ E^ - e , pf ^ pt - k + q . (1) 
The interaction length is determined by gy, Imt ~ 1/9||- From the set of equations 

q±Kp^9, (?|| -cos6l) w ^^,6*2 and 9 ~ hnt/Rc ~ l/{q\\Rc) , (2) 



we find g|| and q± as 



g_L ~ m ( 3 ' ) and g,, « , (3) 
\m-^Rc J 2 Pi 
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where m is the proton mass. 

Naively, the natural parameter of the problem is l/(rni?c)- However, we shall choose as a parameter the quantity 
proportional to this value, namely the quantity in the brackets of Eq. (^): 

X = Ef/{m'Rc) . (4) 

The motivation for this choice is given by the connection with the fictitious magnetic field method. Let us assume 
that the considered part of the trajectory is caused by a fictitious magnetic field B perpendicular to the trajectory 
plane. This means that the Larmor radius Rl = E/{eB) equals the local curvature radius Rc- Then Eq. (^) gives 

X^{E,/m){B^/B,r) . (5) 
where Bcr = vr? je and x is the standard parameter 



X = e^{F^^^lvr? (6) 

relevant for processes in strong electromagnetic fields. In this paper we shall confine ourselves to the quantum case 
X ^ 1 when the curvature radius is small, Rq <C E"^ /m?. 
Two remarks are in order. 

The critical magnetic field in our problem is connected with the proton mass and is equal to B^r ~ 1-5 • 10^° G. 

One should distinguish between the parameter x for the curvature radiation given by Eq. (||) and the parameter 
Xem relevant for the real magnetic field in the considered system. The latter is given by Eq. ^) and results in Eq. (||) 
with i3_L being the component of the real magnetic field B perpendicular to the proton momentum. The difference 
between these two parameters, x given by Eq. (^) and Xem given by Eq. is very essential for a proton moving 
along a curved magnetic field line. The synchrotron radiation of the proton in the real magnetic field results in the 
loss of transverse momentum in the regime when Xem ^ 1^ 



I dp±^ 2 e^m ^ 
p± at 37 



(7) 



where 7 is the Lorentz factor of the proton. For Rc < 10^ cm the lifetime given by ^ in the regime Xem <C 1 is 
shorter than the characteristic time of flight and therefore a proton moves along the magnetic field line At the 
same time the curvature radiation parameter (|^) can be larger (or much larger) than I, providing thus pion production 
due to curvature radiation. In other words, this movement is characterized by Xem 1 and x > 1 (or X S> 1). 

III. SEMICLASSICAL SOLUTION 

Let us consider a proton moving along a curved trajectory. As it is described in the Introduction we single out the 
interaction region according to Eq. (|^) and take the wavefunctions of the proton at the boundaries as plane waves ipi 
and ijjf. Then the matrix element describing the interaction is 

M^g{Pf{pf)^''^,{p,)m . (8) 
where Pi,pf and k are determined by Eq. (|l|). The probability of the transition is calculated in the standard way as 



(2^)^ 1 1^,^,2 dl'k Spf 
2E, 2^ ' {2nf2e {2nf2Ef 



= ^ 77 TKZ^^ 7K-£^_ m -Ef- e)5{p, -p,-k + q). (9) 



In the evaluation of the spinor traces we use the kincmatical relations given by Eqs. (1) and (|). As a result we obtain 

2 _ _1 ^ 2/3 



91 = TTF X'" , (10) 



MnEi 16 E^ 

where as — g^/47r and x is given by Eq. (|^) or in terms of the fictitious magnetic field by Eq. (||). 
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IV. FICTITIOUS MAGNETIC FIELD METHOD 



In this approach we assume that the curved trajectory is caused by a fictitious magnetic field B which is chosen so 
that Rl = Rc- Therefore we now have to calculate pion radiation of a proton in a strong magnetic field. This process 
was already considered in Ref. [|| using the crossed field method |^ (see also Ref. [^). Here, we shall calculate the 
probability of this process describing it as a transition of a proton from the Landau level Ni with polarization = ±1 
into the Landau level Nf due to the interaction Lint = gi^fl^^i't'- The transition probability for a proton from the 
Landau level Ni to N f is given by 

Ti^'^'-^E/^^^^in^fcl-^l'' (11) 

where is the angle between the pion momentum k and the magnetic field. The matrix element squared \M.\'^ is 
given in terms of /-functions. For the very large Landau quantum numbers considered here, we can use an asymptotic 
expansion in terms of modified Bessel functions (for the notations used here and hereafter see Ref. Q): 

In,.n. {x) = ^ ^"-^ VJ^K,/,{z) (12) 



where 



X 



2eB 



(13) 



^ 3 B E^ 

2 Brr m 2 



2^,£/i3/2 

3eB(l - e/Ei) 



(14) 



M=l 



2N,eB 



and 



^J■o = 1 



2N,eB 



(15) 



Furthermore, since pions are emitted mostly in the forward direction, i.e. the transition probability is strongly 
peaked at 6* = 7r/2, we can approximate all functions except functions of 1 — x/xq by their values a,t 9 — Tr/2. Then 
we take the ultrarelativistic limit of the normalization constants of the spinors and obtain 



37r2 



1 



f [(1 - nrf)D' -nil- nTf) D{E - l)]^,^/'K^/,{z)Ky,{z) 



(16) 



[(1 - T,T/)(1 -2E + D^ + E^) - 2n (1 - r.r/) D{E - 1)] ^?Kl,^{z) 



with D ~ m/(2Ef) — m/{2Ei) and E — m/ {AEiEf). Note that because of the pseudoscalar nature of the interaction 
only spin-flip transitions occur in the ultrarelativistic case. 

Using now some identities for the Bessel functions and summing over the polarization states of the proton we obtain 
for the spin averaged emission rate 



3\/37r 



dxK^/^(x) +K2/z(v) 



(17) 



For the case ^ ^ 1 we can use for the functions Ki^i{y) their asymptotic expression K^{y) ~ 2^ ^T{^)y ^ for y — > 0. 
Then, summing over all final Landau levels yields the total decay width F 



22/3r2(|)r(i) ^ 3^^ ^ 

18V3^ ' E,^ ^12.6"^i?, ^ 



2/3 



(18) 
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Our width ^ differs only by a numerical factor from the result of Ref. ||. This factor, 9r(2/3)r(l/3)/(2V37r) w 3.0, 
is however rather large. 

On the other hand, the agreement between our semiclassical solution (|l^) and the Eq. (|8|) is surprisingly good. Both 
methods give not only the same functional dependence on Ei and x, but differ only by a numerical factor 1.8. Because 
of the unccrtainities in the kincmatical relations used in the semiclassical calculation we estimate the uncertainity in 
the decay width by a factor 4—5. However, for two other processes in strong magnetic fields, ^ + + and 
7 — > e+ + e^, we obtained excellent agreement between both methods. 



V. ASTROPHYSICAL APPLICATIONS 



The pion curvature radiation needs protons of extremely high energies. From condition x ^ 1 one obtains from (7) 

E, > 6.5 • 10^i?y^ GeV, 

where Rq — cm is the curvature radius in units of 10^ cm. In cosmic rays the particles are observed with 

energies up to 2 • 10^^ GeV. The particles with the highest energies are most probably protons. In the strong magnetic 
field inside a source protons of these energies move along the magnetic field lines Since the proton momentum is 
essentially parallel to the magnetic field the electromagnetic parameter, Xem ^ 1 and pion production in the strong 
magnetic field is absent. However, the curvature parameter (7) can be large if the curvature radius of the magnetic 
field line is sufficiently small and therefore the pion curvature radiation is efficient. Decays of pions result in the 
production of gamma and neutrino radiation. 
An example can be given by young pulsar. 

Independently from the concrete model of pulsar magnetospheres, the maximum potential 4>max on the surface of 
of the pulsar is 1^: 

^raax > 8.3 • lO^Cfi/lO^ s'^iBjlO^ G)(i?,/106 cm)^ V, (19) 

where Rs is the star radius, il is the angular velocity and Bs is the surface magnetic field. 

Therefore, pulsars with > 3 • 10^ s~^ can, in principle, accelerate protons to energies higher than 1 • 10^° GeV. 
If this potential drop occurs nearby the surface of the pulsar (which is actually not the case of the modern models) 
the ultrahigh energy protons escape from the pulsar magnetosphere mostly along the first-open magnetic field line 

. The curvature radius of this line is 

Rc = 4/3 {RsRLt'^ , (20) 

where Rl = c/Vl is the radius of the light cylinder. For = 3 • 10'^ s~^ and Rs = 1.2 • 10^ cm, one obtains from Eq. 
( p^ that the energy of accelerated protons could be higher than 1.3 • 10^" GeV and the parameter (^ is x > 1. 

Another example is an accretion disc, with a strong magnetic field, around a black hole with a mass M ~ M©. The 
potential at the inner edge of the disc is |Q 

1 /" 1 /" 1 R 

(j)=- Edl = - / v{r)B^{r)dr = — B^ro In - , (21) 
c J c J V6 rQ 

where rg, R and Bq are the radii of the inner and outer edges of the disc and the magnetic field at the inner edge, 
respectively. Taking Ei ~ e0, tq ^ Rc ^ 10^ cm and Bq ^ 10^° G we again obtain x > 1. 

Finally, let us discuss superconducting cosmic strings. Cosmic strings are topological defects caused by symmetry 
breaking in the early universe. Under certain assumptions the string can become superconducting Due to the 
expansion of the universe the primordial magnetic field induces an electric current J in the string which grows up to 
the maximal value Jmax ~ eM/27r, where M is the scale of the GUT symmetry breaking. Starting from this moment, 
the particles with GUT masses are evaporating from the string. Decays of these particles produce leptons and hadrons 
with maximum energies up to Mqut ~ 10^^ GeV. 

The string current produces a magnetic field with Biot-Savart geometry. The initially produced particles have 
momenta perpendicular to the magnetic field. Therefore, at a first stage the processes in the strong magnetic field 
dominate . After losing energy, particles are captured in circular orbits |lO[| and the pion curvature radiation becomes 
one of the main energy loss mechanisms of high energy protons. 
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VI. CONCLUSIONS 



We have studied the curvature radius of pion by a proton. The parameter of this problem is x = Ef / (m^Rc), where 
Rc is the curvature radius of the proton trajectory and Ei is the proton energy. 

For the case x ^ 1 we found the quasiclassical solution assuming that in the region of coherent interaction the 
trajectory gives the additional momentum q to the proton. The parallel component of this momentum, gj| , determines 
the length of coherent interaction lint ~ Then the perpendicular momentum, g_L, is found as g_L ~ PiO, where 

Pi is the proton momentum and 9 ~ hnt/Rc- The wavefunctions of the proton at the boundaries of the interaction 
region, Unt, are taken as plane waves and the matrix element is calculated as gu{pf)"fQu{pi)(j) with pf and pi related 
as Pf ^ Pi — k + q, where k is the pion momentum. The width of the proton decay, p ^ p + ir'^ is given by 

r « 0.19asm2x^/V-E^ . (22) 

We estimate the uncertainities of this method by a factor 4 — 5. 

The other way of calculation is given by the method of fictitious magnetic field. Here we assume that the proton 
moves along the trajectory, due to a magnetic field, so that locally Rl = Rc, where Rl is the Larmor radius of 
the trajectory. We calculated the width of proton decay in this case as a probabilty of a transition between different 
Landau levels and found 

r w 0.104 a^TO^ ^2/3/^^ (23) 

in reasonable agreement with the quasiclassical method. 

However, the early calculations of Rcf. Q for the p ^ p + tt" decay in the strong magnetic field, made using a 
different approach, differ from our calculations by factor 3. 

As applications we consider the astrophysical sources where high energy protons are moving along magnetic field 
lines. Since the proton momentum is essentially parallel to the magnetic field, the electromagnetic parameter Xem 1, 
while the curvature radiation parameter can be x > 1 (in the case of sufficiently small Rc)- These conditions are met 
in such sources as young pulars, accretion dies around black holes, and superconductive cosmic strings. 
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